Anisotropic flow of charged and identified hadrons in the quark-gluon string model 

for Au+Au collisions at ^/snn = 200 GeV 
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The pseudorapidity behaviour of the azimuthal anisotropy parameters vi and V2 of inclusive 
charged (h ± ) hadrons and their dependence on the centrality has been studied in Au+Au collisions 
at full RHIC energy of ^snn = 200 GeV within the microscopic quark-gluon string model. The 
QGSM simulation results for the directed flow V\ show antiflow alignment within the pseudorapidity 
range \n\ < 2 in a fair agreement with the experimental Vi(rf) data, but cannot reproduce the further 
development of the antiflow up to \r]\ ~ 3.5. The n dependence of the elliptic flow V2 extracted from 
the simulations agrees well with the experimental data in the whole pseudorapidity range for different 
centrality classes. The centrality dependence of the integrated elliptic flow of charged hadrons in the 
QGSM almost coincides with the PHOBOS experimental distribution. The transverse momentum 
dependence of the elliptic flow of identified (tt , K , p, p) and inclusive charged hadrons is studied 
also. The model reproduces quantitatively the low-pr part of the distributions rather good, but 
underestimates the measured elliptic flow for transverse momenta pr > 1 GeV/c. Qualitatively, 
however, the model is able to reproduce the saturation of the V2(pr) spectra with rising pr as well 
as the crossing of the elliptic flow for mesons and baryons. 

PACS numbers: 25.75.-q, 25.75.Ld, 24.10.Lx, 12.40.Nn 

Keywords: ultra-relativistic heavy-ion collisions; directed and elliptic flow; pseudorapidity and transverse 
momentum dependence; Monte-Carlo quark-gluon string model 
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I. INTRODUCTION 

Ultra-relativistic heavy ion collisions (URHICs) offer a 
unique opportunity to study the nuclear phase diagram 
at high temperatures and densities 0. The matter un- 
der such extreme conditions has probably existed in the 
early universe within the first few fm/c's after the Big 
Bang. Therefore, it is very tempting to investigate the 
properties of the Little Big Bang [2j in the laboratory, 
and to search for a new state of matter, predicted by 
the fundamental theory of strong interactions (Quantum 
Chromodynamics - QCD), namely, a plasma of decon- 
fined quarks and gluons (QGP). 

Among the various experimental studies of URHICs 
operates the Relativistic Heavy Ion Collider (RHIC) at 
BNL since 2000 to investigate gold on gold collisions 
up to y/sNN = 200 GeV. After four years of operation 
strong experimental evidence has been accumulated that 
at RHIC energies indeed a new state of matter is created 
which is qualitatively different from a hadron gas (see Q 
and references therein). This state seems, however, not 
to behave like a weakly interacting parton gas - as could 
have been naively expected - but rather like a strongly 
coupled plasma. One argument towards such a scenario is 
the large elliptic flow observed at RHIC 0, H H, H . The 
development of strong elliptic flow requires short equi- 
libration times and large pressure gradients to drive the 
dynamics, both being characteristic features of a strongly 
interacting system. 

An independent argument for such a scenario is pro- 
vided by lattice QCD, though still on a qualitative ba- 



sis. While lattice QCD predicts undoubtedly a QGP 
phase transition around a critical temperature of T c ~ 
150 -=- 180 MeV the properties of such a state are not 
yet so well understood. As a striking fact lattice calcu- 
lations do not reach the Stefan-Boltzmann limit for the 
pressure of a free parton gas p oc T 4 even at T > 5T C but 
saturate around 80% of the Stefan-Boltzmann limit • 
Taking this fact seriously and not as a lattice artefact, 
one may argue that the QGP at temperatures well above 
the critical one is still a strongly interacting system || . 

In line with these arguments goes the success of hydro- 
dynamics in describing the collective flow at RHIC. Only 
at incident energies of ^/snn > 130 GeV (which corre- 
sponds to an initial energy deposit of about 10 GeV/fm 3 
in the central reaction zone) the hydrodynamic limit for 
the elliptic flow seems to be reached . Within the limits 
of tuning of the equation of state (including incorpora- 
tion of a QGP phase transition) hydrodynamic calcula- 
tions are able to describe the bulk properties of the col- 
lective flow rather reasonable [lflllllll2fl . However, when 
turning to more differential observables such as centrality 
dependence, mass dependence or px dependence of the 
elliptic flow, hydrodynamic calculations have also some 
problems to match the data. 

In this context it is also important to obtain an un- 
derstanding of the reaction dynamics in terms of micro- 
scopic models. One class of microscopic models which 
have been very successfully applied at CERN SPS ener- 
gies and below are string models. The main assumption 
of string models is that hadrons are produced as a re- 
sult of excitation and decay of open strings with different 
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quarks or diquarks at their ends. Generally, all models 
are formulated as Monte Carlo event generators allowing 
to perform a careful analysis of the measurable quantities 
by introducing all necessary experimental cuts. There 
are numerous versions of the two basic string-motivated 
phenomenological approaches: the FRITIOF model 0] 
and the Dual Parton Model (DPM) 14]. These two ap- 
proaches use different mechanisms of string excitation. 
In the first approach, relied on relativistic classical string 
dynamics, the string masses arise from momentum trans- 
fer. The second one is based on the Gribov-Regge the- 
ory (GRT) 15] in the framework of relativistic quantum 
theory where quantum aspects like unitarity play an es- 
sential role. Here, the strings are produced as a result of 
colour exchange. 

By construction, such types of models do not con- 
tain explicitly a quark-hadron phase transition. However, 
during the temporal evolution of a heavy ion reaction a 
dense and strongly interacting plasma is formed within 
such approaches as well. The system consists of partons 
and colour- flux tubes (or strings). Thus, it is an essen- 
tial question if string models are able to create a suffi- 
cient amount of pressure in order to produce large elliptic 
flow seen at RHIC, which features they can describe and 
where they might fail. Such investigations are in particu- 
lar relevant since transverse as well as elliptic flow at SPS 
energies is well rep roduced within the string-cascade ap- 
proach [l(| [T3, d, EH • By such type of studies one can 
obtain deeper insight in the question which observables 
indicate the appearance of new physics not included in 
standard approaches to relativistic heavy ion collisions. 

In the present work we describe ultra-relativistic heavy 
ion reactions by a microscopic quark-gluon string cascade 
model (QGSM) based on Gribov-Regge theory. Details 
of the model are given in the next section. The QGSM 
has been demonstrated [2(| to give a fair description of 
the first data of v-i in Au+Au reactions at ^/snn = 130 
GeV kJ. Here, we extend the analysis to full RHIC en- 
ergy (^snn — 200 GeV) and to more differential observ- 
ables, i.e. centrality dependence, mass dependence and 
Pt dependence of V\ and V2- 



II. QUARK-GLUON STRING MODEL 

The Quark-Gluon String model (QGSM) incorporates 
partonic and hadronic degrees of freedom and is based on 
GRT accomplished by a string phenomenology of parti- 
cle production in inelastic hadron-hadron (hh) collisions. 
To describe hadron-hadron, hadron-nucleus and nucleus- 
nucleus collisions the cascade procedure of multiple sec- 
ondary interactions of hadrons was implemented. The 
QGSM incorporates the string fragmentation, formation 
of resonances, and rescattering of hadrons, but simpli- 
fies the nuclear effects neglecting, e.g., the mean fields or 
evaporation from spectators. As independent degrees of 
freedom the QGSM includes octet and decuplet baryons, 
octet and nonet vector and pseudoscalar mesons, and 



their antiparticles. The momenta and positions of nu- 
cleons inside the nuclei are generated in accordance with 
the Fermi momentum distribution and the Woods-Saxon 
density distribution, respectively. Pauli blocking of oc- 
cupied final states is taken into account. Strings in the 
QGSM can be produced as a result of the colour exchange 
mechanism or, like in diffractive scattering, due to mo- 
mentum transfer. The Pomeron, which is a pole with 
an intercept ap(0) > 1 in the GRT, corresponds to the 
cylinder-type diagrams. The s-channel discontinuities of 
the diagrams, representing the exchange by n-Pomerons, 
are related to the process of 2k (k < n) string produc- 
tion. If the contributions of all n-Pomeron exchanges 
to the forward elastic scattering amplitude are known, 
the Abramovskii-Gribov-Kancheli (AGK) cutting rules 
[2l| enable one to determine the cross sections for 2k- 
strings. Hard gluon-gluon scattering and semi-hard pro- 
cesses with quark and gluon interactions are also incor- 
porated in the model [22j. The inclusive spectra in the 
QGSM have automatically the correct triple- Regge limit 
for the Feynman variable x — ► 1, double- Regge limit for 
x — > 0, and satisfy all conservation laws |2^]. The partic- 
ular stages of the collision model, namely (i) initializa- 
tion of interacting projectile and target nuclei, (ii) string 
formation via inelastic nucleon-nucleon (hadron-hadron) 
interaction, (hi) string fragmentation, i.e. hadronization, 
and (iv) hadron-hadron rescattering, are solved basically 
by Monte Carlo simulation techniques |24j . 



III. ANISOTROPIC FLOW OF INCLUSIVE 
CHARGED AND IDENTIFIED HADRONS 

Among the main signals, which can help to reveal the 
formation of the QGP in the experiment, are collective 
flow phenomena. Flow is directly linked to the equation 
of state of the excited matter. Generally, an effective 
EOS extracted from the model studies shows "softness" 
during the early stages of the collision, but it remains 
unclear whether the observed softness is (i) due to the 
proximity of the QCD phase transition |25l l26l 1271 128| . 
or (ii) due to non-equilibrium phenomena, such as the 
formation and fragmentation of strings 

jMlSJ, or (in) 

due to nuclear shadowing [Til WK Isil ] . 

The transverse collective flow can be subdivided into 
isotropic and anisotropic flow. Two types of anisotropic 
flow, that are characterized by the first and the second 
harmonic coefficients of the Fourier decomposed invariant 
azimuthal distribution in momentum space 32, 33] 
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are called directed and elliptic flow. Here, pt = (j> 2 + 
Py) 1 / 2 is the transverse momentum, y the rapidity and 
<p the azimuthal angle between the particle's momen- 
tum and the reaction plane. While the elliptic flow 
v 2 = (cos(20)) = ((px/pt) 2 — (Py/pr) 2 ) is assumed to be 
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particularly sensitive to the pressure at the early stages 
of the collisions [l(J, 123, H3, HH , the directed flow appears 
to develop until the late stage of the heavy-ion reaction 
\wi ll7L H8| . On the other hand, the directed flow of par- 
ticles with high transverse momentum, which are emitted 
at the onset of the collective expansion, can carry infor- 
mation about the EOS of the dense matter phase from 
the initial conditions. A lot of measurements of the vi 
parameter in Au+Au collisions at RHIC have beenper- 
formed for charged and identified hadrons, see [HOHEE! 
and references therein, whereas experimental results for 
the directed flow v\ at RHIC has been reported only quite 
recently plE^]. 



A. Directed flow of inclusive charged hadrons 

Figure depicts the QGSM simulation result for the 
pseudorapidity dependence of the directed flow V\(j]) of 
charged hadrons compared to the experimental data from 
the PHOBOS Collaboration 36] for 6% to 55% central 
Au+Au collisions at ^JJWn = 200 GeV (QGSM pre- 
dictions for the directed flow of identified particles at 
lower RHIC energy, ^snn = 130 GeV, can be found in 
[37|]). One can see that the model reproduces the V\ 
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FIG. 1: Directed flow vi for charged hadrons as a function of 
pseudorapidity r\ in comparison to the result from the PHO- 
BOS Collaboration (centrality 6% to 55%) 36] for Au+Au 
collisions at ^/snn = 200 GeV. Statistical (bars) and sys- 
tematic (boxes) experimental errors are shown. 



tic wiggle structure with a clear antiflow 1 component in 
the middle |?7| region. It was pointed out in |3^ (see 
also 0, Hil) that the phenomenon leading to the for- 
mation of a wiggle structure for the directed flow of nu- 
cleons is caused by dense baryon-rich matter shadowing, 
which plays a decisive role in the competition between 
normal flow and antiflow in noncentral nuclear collisions 
at ultra-relativistic energies. Within microscopic string 
model calculations such deviations from the straight line 
behaviour of the nucleonic flow were first observed in 
very peripheral Au+Au collisions at AGS energy [3l| . 
Experimentally the wiggle structure of v\ for protons in 
peripheral Pb+Pb collisions at SPS has been observed by 
the NA49 Collaboration USE!. However, the QGSM 
distributions of vi(rj) at these lower energies also have 
peaks which are shifted b y appr oximately one unit of 
rapidity toward 77 = [l6l Il7l 12^] compared to the ex- 
perimental data. So, this shift, which also leads to a 
steeper slope in the mid-rapidity region, seems to be 
a sort of model 'artefact' occurring at higher collision 
energies as well. The question about the origin of the 
directed flow's rapidity dependence obtained within the 
QGSM is not so easily to answer. One reason might 
be the lack of heavy resonances in the model compared, 
e.g., to the Ultra-Rclativistic Quantum Molecular Dy- 
namics (UrQMD) model ^3] or the multi-phase trans- 
port (AMPT) model 0]. These resonances are mostly 
concentrated in the rapidity regions of flying away specta- 
tors which exhibit normal flow. Lighter particles coming 
from the decays of the heavy resonances, each moving in 
the normal direction and having nearly the same rapid- 
ity, may significantly enhance the normal component of 
the directed flow in a pseudorapidity range closer to the 
fragmentation region. The total multiplicity of particles 
with pseudorapidity \rj\ > 3.5 is quite low in the present 
version of the QGSM. The complex connection between 
the model dynamics of the QGSM and the characteristics 
of the resulting directed flow requires further lucid inves- 
tigations to understand in detail the origin of the change 
in sign, and in particular the strong antiflow behaviour 
of V\ (77) in the mid-rapidity region. 

Nevertheless, our picture is a bit different compared 
to that provided by the microscopic models based on the 
FRITIOF routine. For instance, the relativistic quantum 
molecular dynamics (RQMD) model favours weak but 
still normal flow for pions |38| even for more peripheral 
topologies with b = 5 — 10 fm, corresponding to a central- 
ity range a/a geo = 15 — 60%. Both the UrQMD model 
and the AMPT model show a very flat and essentially 
zero directed flow 0, 03 m a broad range \r]\ < 2.5. 
- It is worth mentioning here that the results of both 
models have been obtained for minimum bias, not semi- 



data quite well both qualitatively and quantitatively, al- 
though the maxima of the directed flow around \r)\ 2 
are shifted to lower pseudorapidities compared to the ex- 
perimental data. The vi{rf) result show a characteris- 



Conventionally, the type of flow with positive slope dvi/dr) is 
called normal flow, in contrast to antiflow for which dvi/dn is 
negative. 
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FIG. 2: Directed flow v\ for nucleons and charged pions from y^ijvjv = 200 GeV Au+Au reactions as a function of pseudorapidity 
n for six different impact parameters b going from central (6 = 2 fm) to peripheral (b = 12 fm) collisions. Only statistical errors 
are shown. 



peripheral, events. - Although the data seem to indicate 
antiflow behaviour for the directed flow of charged parti- 
cles with the possible flatness at |f?| < 1.5, the measured 
signal is quite weak, - the magnitude of the flow is less 
than 1% at \r)\ < 2. Therefore, relatively large system- 
atic error bars do not permit to disentangle between the 
different models. 

Similar antiflow alignment can be obtained also within 
the multi module model (MMM) ,441] , which is based on 
fluid dynamics coupled to the formation of colour ropes. 
In this model the effect of the tilted initial state, respon- 
sible for the antiflow formation, reaches its maximum for 
the impact parameter b ~ 0.5(Ra + Rb), i-e. for the 
centrality o~/o~ geo ~ 25% in case of a symmetric system 
of colliding nuclei. To check the centrality dependence 
of the directed flow at full RfflC energy, the pseudo- 
rapidity distributions vi(n) of nucleons and charged pi- 
ons in Au+Au collisions at ^snn = 200 GeV are dis- 
played in Fig. |21 for six different impact parameters b 
going from central (upper left panel) to very peripheral 
(lower right panel) configurations. Regardless of the cen- 
trality range, directed flow of pions has negative slope, 
i.e. antiflow, in the mid-rapidity range. For nucleons the 
azimuthal anisotropy parameter v\ is small for central 
collisions (6 = 2 fm). In accord with our previous stud- 
ies and conclusions ll6l for Au+Au collisions at full 
RHIC energy deviations of nucleonic flow from a normal 
flow behaviour occur already at quite small impact pa- 
rameters. This means that the effect is indeed shifted 
to more central configurations. With the increase of b 
these deviations, representing the wiggle structure of the 
flow, appear more distinctly. However, the QGSM simu- 
lations for both nucleonic and pionic flow at mid-rapidity 
\r)\ < 1 are consistent with the zero flow signal v± = 0. 
The other features which should be mentioned here are 



the broadening of the antiflow region and the increase 
of its strength as the reaction becomes more peripheral. 
Therefore, one can disentangle between two processes of 
different origin employed for the description of the third- 
flow component: If the formation of nucleonic antiflow is 
dominated by the creation of QGP, the flow maximum 
is reached at b < 6 fm 44], whereas for the shadowing 
scenario the strong antiflow should be observed also in 
very peripheral events with b ps 10 fm. 



B. Elliptic flow of charged and identified hadrons 

Here we investigate the pseudorapidity dependence of 
the elliptic flow of charged hadrons. The QGSM sim- 
ulation result for minimum bias Au+Au collisions at 
y?SNN — 200 GeV are compared in Fig. |31 with the ex- 
perimental data of the PHOBOS Collaboration [5]. The 
elliptic flow displays a strong in-plane alignment in accor- 
dance with the predictions of Ref . |3J] . At mid-rapidity 
\n\ < 1.0 the elliptic flow is almost constant. Then it rises 
up slightly and drops rapidly at | -^7 1 > 2.0 with increas- 
ing pseudorapidity. The mean value of V2 h (M < 2.5) is 
practically as large as the value measured by the PHO- 
BOS Collaboration at mid-rapidity. The experimental 
data, which indicate a steady decrease in vi with increas- 
ing 1 77 1, are slightly overestimated by the model only at 
1 77 1 ps 2.0 as a consequence of a double hump structure in 
the theoretical result. This difference in shape close to 
r\ = still rankles somewhat, although within the error 
bars the pseudorapidity dependence of the elliptic flow 
of charged hadrons obtained within the QGSM shows a 
really fair agreement with the PHOBOS data p| in the 
whole r\ range. 

However, following the idea of longitudinal boost in- 
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FIG. 3: Elliptic flow «2 for inclusive charged hadrons as a 
function of pseudorapidity r\ in comparison to the PHOBOS 
data of minimum bias Au+Au collisions at y^ijvjv = 200 GeV 
0- The systematic errors of the experimental data are shown 
as gray boxes together with the statistical errors (bars), re- 
spectively. 



variance 2 of the expanding hot and dense matter and the 
common interpretation of elliptic flow as a consequence 
of secondary particle collisions, one would expect no or at 
least a weak pseudorapidity dependence of v-i similar to 
that of the multiplicity density dN / dr\ over a large rapid- 
ity range. The experimentally observed multiplicity stays 
approximately constant within three units of pseudora- 
pidity |4i| , while the elliptic flow data show a pronounced 
peak at mid-rapidity [H ^(J . This is somehow in contra- 
diction with the assumption of longitudinal boost invari- 
ance over a broad region of rapidity in RHIC collisions. 
Most of the hydrodynamics calculations reported in the 
literature are based on boost invariant models. Therefore 
the results obtained within such approaches are indepen- 
dent of rapidity and one is limited to discuss only the 
transverse behaviour • Hydrodynamics results for the 
pseudorapidity dependence of v 2 are scarce. The shape 
of 1*2(77) in Au+Au collisions at ^/snn = 130 GeV stud- 
ied within a full three-dimensional hydrodynamic model 
by Hirano et al. [iM E^ . IHof shows a bump structure in 
forward and backward rapidity regions, when the initial 
energy density profile depends strongly on the space-time 
rapidity r/ s . The bumps appear at the same rapidity re- 
gions where the initial configuration is not boost invari- 
ant anymore. In the case of an almost r/ s independent de- 
formation of the initial energy density, the resultant i>2 (rj) 



2 The assumption of a longitudinal boost invariant system means 
that its energy density and pressure do not depend on the 
longitudinal coordinate z compared at the same proper time 
r = Vi 2 — z 2 . In other words, the evolution of the pressure 
and energy density depends only on r but not on r\. 



has no bumps. Thus it was shown in [50j that the pseudo- 
rapidity dependence of the elliptic flow is highly sensitive 
to the parametrization of the initial energy density profile 
in the longitudinal direction. The initial energy density 
profile in microscopic transport models is not explicitly 
parametrized, but it is implicitly fixed by the initial con- 
ditions of the projectile and target nucleus. The elliptic 
flow studied within the UrQMD model in the cascade 
mode shows also a prominent dip at central rapidities for 
all inspected hadrons |4^. There it is argued that this 
rapidity behaviour of V2 indicates a region of small inter- 
action strength (or low 'pressure') because of the direct 
connection between the strength of the anisotropic flow 
and the mean free path of the particles forming the hot 
mid-rapidity region. Therefore, the appearance of the 
dip in V2 is linked to a feature of the model dynamics 
in the early stage, namely the pre-equilibrium string dy- 
namics and interactions on the parton level |4^. Also in 
the QGSM, the emergence of the double bump structure 
in 172(17) is strongly connected with the model dynamics, 
as one can clearly see in the time evolution of the elliptic 
flow and its rapidity dependence |5l|. The normalized 
number of binary particle collisions per pseudorapidity 
interval throughout the evolution of the system as a func- 
tion of the pseudorapidity is depicted in Fig. 0] for two 
impact parameters b — 2 fm and b = 8 fm to show the 
strong correlation between the particle interactions and 
the flow strength mainly produced by secondary colli- 
sions. First of all, this number versus the pseudorapidity 
difference A77 = 771 — 772 of the two colliding particles 
is plotted in the left panel of Fig. 0] The peak around 
A?; = indicates unambiguously that independently of 
the impact parameter the particles with equal or at least 
very similar rapidities interact most likely. The plot on 
the right hand side of Fig. 01 where the number of binary 
collisions normalized to its maximum is shown as a func- 
tion of the mean pseudorapidity fj = (771 + 772) /2 of the 
two colliding particles, is even more instructive. Here, a 
clear double peak structure appears in the same rapidity 
region as seen in Fig. [3] for the 77 dependence of the ellip- 
tic flow. The dip in the number of collisions around fj = 
becomes more pronounced with decreasing centrality. In 
addition, the two peaks for b — 8 fm are slightly shifted 
to higher pseudorapidities compared with the result for 
6 = 2 fm. This is nicely consistent with the centrality de- 
pendence of V2 (??) obtained within the QGSM as it will 
be discussed in the following. 

Figure [3] presents the pseudorapidity distribution of 
the elliptic flow for charged hadrons in gold-gold colli- 
sions at full RHIC energy for three different centrality 
classes, ranging from central (bottom panel) via mid- 
central (middle panel) to peripheral (top panel) in ac- 
cordance with the definitions in [5j . The results from the 
QGSM simulation and the PHOBOS analysis (combined 
data from the hit- and track-based methods) are overlaid. 
The model is able to describe the magnitude and shape 
of vf 1 ^) quite well across all of the three centrality bins 
within the given systematic and statistical uncertainties. 
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FIG. 4: The number of binary particle collisions per pseudorapidity interval normalized to its maximum as function of the 
pseudorapidity difference A77 = 771—772 (left) and the mean pseudorapidity 77 = (771 + 772)/2 (right) of the two colliding particles. 
The QGSM results are obtained for simulated Au+Au reactions at -^/sjvjv = 200 GeV with two different impact parameters b. 
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FIG. 5: Elliptic flow ?J2 for inclusive charged hadrons from 
^/sjvjv = 200 GeV Au+Au collisions as a function of pseu- 
dorapidity 77 for the three centrality classes according to the 
PHOBOS analysis (combined hit- and track-based results) 
[f|- Again, the systematic errors of the experimental data are 
shown as gray boxes together with the statistical error bars, 
respectively. 



The overall shape of the 77 distribution changes only very 
little with centrality and shows a behaviour very similar 
to that depicted in Fig. [31 for the minimum bias events. 
The two-peak structure of vf 1 ^) in the interval |?7| < 2.5, 
which is clearly seen for simulated peripheral collisions 



but 'washed out' in the central bin, is not attributed 
solely to QGSM, but arises also in e.g. UrQMD calcu- 
lations |42| as discussed in detail above. However, the 
effect is small, and the mean value of the vf 1 parameter 
over the aforementioned range increases from central to 
peripheral collisions in good quantitative agreement with 
the experimental data. This is a not so trivial result, 
because neither the pseudorapidity nor the centrality de- 
pendence of the elliptic flow, which is discussed below, 
is reproduced correctly at RHIC so far by the UrQMD 
calculations and hydrodynamic models. 

The next observable is the transverse momentum de- 
pendence of V2 for identified hadrons, namely, combined 
7r + + 7r~ , K + + K~ , as well as p + p spectra in minimum 
bias Au+Au collisions at ^/sWn = 200 GeV. The QGSM 
simulation results are depicted in Figs.|Hland0in com- 
parison to the PHENIX |(J and STAR data [Bj, respec- 
tively. The agreement with the experimental data at least 
for protons/antiprotons and kaons is rather good in the 
range of low transverse momenta 0.0 < pr < 1.0 GeV/c, 
as illustrated in Fig.El The steady increase of the elliptic 
flow with rising px and the larger vi parameter at a given 
Pt for the lighter mass hadrons compared to the heavier 
ones is well reproduced by the model calculation. For 
pions the latter statement is only valid for the transverse 
momenta below 0.75 GeV/c. The elliptic flow of charged 
pions starts already to saturate at pr > 0.5 GeV/c. This 
early saturation of the pionic vi continues with a slight 
decreasing at increasing px and entails the peculiar result 
to be smaller than the kaon flow for transverse momenta 
between 1.0 GeV/c and 2.5 GeV/c, as shown for the over- 
all pt range in Fig. [7] 

This deviation compared to the experimental findings 
is remarkable because it was shown in 20] that the mag- 
nitude of the pionic flow in the QGSM calculations is 
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FIG. 7: The same like in Figure HJ but in the overall pr range 
0.0 < p T < 3.5 GeV/c. The PHENIX data are taken from 
0- Again, only statistical errors are shown. 



already twice as large as obtained, e.g. in the RQMD 
ones. It has been elaborated in [2(| that the contribu- 
tions of hard processes and multi-Pomeron exchanges are 
very important to reach the reported magnitude of v 2 and 
to reproduce the particle multiplicities correctly. To con- 
firm the latter statement, the transverse mass spectra of 
positively and negatively charged pions and kaons for pp 
collisions at ^/snn — 200 GeV are presented in Fig. |S] 
The rriT spectra obtained within the full version of the 
QGSM, which incorporates hard and multichain contri- 
butions, are in good agreement with the experimental 
spectra reported by the STAR Collaboration 
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FIG. 8: Transverse mass spectra of 7T and K for simulated 
pp collisions at ^/snn = 200 GeV compared to the STAR 
data 1531. 



Coming back to the pr behaviour of the elliptic flow 
depicted in Fig-Hit is important to stress at this point an- 
other striking feature of the QGSM. Namely, this model 
is able to reproduce at least qualitatively the experimen- 
tal evidence of the crossing of the elliptic flow for mesons 
and baryons at px ~ 1.7 GeV/c [6j. The simulation data 
show that at pr < 1-4 GeV/c the v 2 parameter for kaons 
is larger than for (anti)protons. At higher transverse 
momenta px > 1-4 GeV/c the situation is completely 
changed. Here, the elliptic flow of protons and antipro- 
tons becomes larger than vi for charged kaons and pions. 

This behaviour is also described within the quark re- 
combination models |54l I5H l56l IS?} , which assume the 
statistical coalescence of two or three quarks into a 
hadron. In contrast, the hydrodynamic picture shows 
the same mass-ordering for the elliptic flow of different 
particles at all transverse momenta [llj . We will come to 
this point after the study of vi(j)t) for inclusive charged 
hadrons. The QGSM distributions are shown in Fig. 
in comparison to the experimental data of the PHOBOS 
and PHENIX Collaborations HQ- To demonstrate the 
importance of rescattering processes for the elliptic flow 
formation the v% h extracted from a QGSM simulation run 
for minimum bias Au+Au collisions without the hadronic 
cascade is additionally plotted in this figure. It is quite 
obvious that the model without subsequent secondary in- 
teractions of produced hadrons creates zero elliptic flow, 
because the azimuthal distributions of secondaries in el- 
ementary hadron-hadron collisions are isotropic. The 
anisotropic flow, i.e. the azimuthal anisotropics of the 
number of produced hadrons, results from the spatial 
asymmetry of the collision zone and subsequent rescat- 
tering processes, which are crucial to convert the initial 
spatial anisotropy into the final momentum anisotropy. 
Hadronic rescattering including hard and multichain con- 
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tributions creates an elliptic flow, which rises almost lin- 
early according to the experimental data within the in- 
terval 0.0 < pt < 1.0 GeV/c, but saturates already 
for transverse momenta above 1.0 GeV/c at a level of 
vf 1 w 6% whereas the experimentally measured flow in- 
creases further up to vf 1 w 14 — 16% ere it declines for 
Pt > 3.0 GeV/c. Although the model seems to describe 
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the pt dependence qualitatively well, it underestimates 
the experimental elliptic flow of charged particles with 
transverse momenta above 1.0 GeV/c roughly by up to 
50%. What is the reason of these deviations? Recall that 
hadrons in the QGSM gain the transverse momentum 
due to (i) transverse motion of the constituent quarks 
and (ii) transverse momentum of constituents acquired 
in the course of string fragmentation. The parameters 
of these two processes are fixed by comparison with the 
available hadronic data. Other sources of the transverse 
motion are (iii) the transverse Fermi motion of nucleons 
in the colliding nuclei and (iv) rescattering of produced 
particles in the hot and dense nuclear medium. It was 
already mentioned that the latter process is the most 
crucial for the development of the elliptic flow in nuclear 
reactions. However, secondary hadrons with high trans- 
verse momenta experience in average less collisions than 
their low-momentum counterparts because of the large 
formation time (which originates from the uncertainty 
principle) . As was shown in |42| , in the limit of vanishing 
formation time the elliptic flow increases drastically. But 
all parameters linked to the formation time of produced 
particles in the QGSM are also fixed by comparison with 
experimental data on hadronic interactions. Therefore, 
the presented elliptic flow can be considered as an upper 
limit obtained within the hadronic cascade scenario. It 
clearly indicates on new physical effects not taken into 



account by the microscopic model. Jets are among the 
most likely candidates for these processes. Indeed, as 
was discussed in [S^ , the non-uniform dependence of the 
energy loss on the azimuthal angle results in azimuthal 
anisotropy of jet spectra in non-central nuclear collisions. 
This leads to a significant increase of the elliptic flow of 
high-pr particles. Another possible explanation is, e.g. a 
dramatic increase of all s-channel transition rates in the 
vicinity of the chiral phase transition |5^|, which causes 
a critical opacity and fast thermalization in the system. 

Finally, the centrality dependence of the total, i.e. rj 
and pt integrated V2 for charged hadrons is depicted in 
the middle and right panel of Figure ^| Since the cen- 
trality of collision events in the experimental measure- 
ments is characterised by the mean number of par- 
ticipants A par t (seen in the middle panel), we show in 
addition to this signal the original impact parameter de- 
pendence of the elliptic flow (right panel). The relation 
between the number of participants A par t in the Au+Au 
reactions and the impact parameter b of the simulated 
collisions is shown in the left panel of Figure ^| The 
number of nucleons participating in inelastic scatterings, 
i.e. Ap art , is directly available in the QGSM. This num- 
ber depends mainly on the nucleon-nucleon cross section, 
which for the given initial energy of ^snn — 200 GeV is 
calculated in the framework of GRT. Apart from the cross 
section and its energy dependence, the number of partici- 
pants varies with the impact parameter like cx exp(— b 2 ). 
Thus, the 6-dependence of A part within the QGSM is 
rather different compared to the Glauber model result, 
but in nice agreement with the experimentally estimated 
Ap ar t (6) H3- Amazingly, the elliptic flow of inclusive 
charged hadrons extracted from our simulation almost 
coincides with the PHOBOS experimental data. This 
remarkable result in conjunction with the transverse mo- 
mentum dependence of vf 1 shown in Fig. reflects the 
dominance of hadrons with low pt- Furthermore one 
can see that the elliptic flow as a function of the impact 
parameter b reaches a maximal value of around 6% at 
b w 8 fm, corresponding to A par t ~ 100, and decreases 
with further increasing b. As expected, the flow in the 
mid-rapidity region is caused mainly by pions, but the 
azimuthal anisotropy parameters Vi (b) for charged kaons 
and combined protons/antiprotons are also very similar 
within the statistical uncertainties. Note that the magni- 
tude of the total pionic flow obtained within the QGSM 
simulations is more than twice as large as created by the 
RQMD ones Q. 



IV. SUMMARY AND CONCLUSIONS 

In summary, the pseudorapidity distributions of the 
azimuthal anisotropy parameters vi(j]) and V2{r)) of in- 
clusive charged hadrons and their centrality dependence 
has been studied in Au+Au collisions at full RHIC en- 
ergy of y/SNN — 200 GeV within the microscopic quark- 
gluon string model. The QGSM simulation results for the 
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directed flow show antiflow alignment within the pseu- 
dorapidity range \r)\ < 2 in a fair agreement with the 
experimental vi(j]) data, but cannot reproduce the fur- 
ther development of the antiflow up to « 3.5. In a 
broad pseudorapidity region the model generates a wig- 
gle structure for the directed flow of nucleons . At 
mid-rapidity < 1, however, the generated flow is quite 
weak and consistent with a zero-flow behaviour reported 
by the STAR and PHOBOS collaborations. The rj de- 
pendence of the elliptic flow v-i extracted from our sim- 
ulation agrees well with the experimental results in the 
whole 1] range for minimum bias as well as for central, 
mid-central and peripheral collisions. The transverse mo- 
mentum dependence of the elliptic flow V2{pt) of iden- 
tified (7r ± , K ± 1 p,p) and inclusive charged hadrons has 
been investigated within the QGSM also. The descrip- 
tion by the quark-gluon string model is fairly good in the 
low pt range. Here, it was shown that for identified and 
charged hadrons the i>2 parameter rises with increasing 
Pt according to the experimental data. For higher trans- 
verse momenta px > 1 GeV/c it starts rapidly to satu- 
rate already on a level, which is at the largest transverse 
momenta roughly 50% smaller than the experimentally 
measured vi- On the other hand, the qualitative be- 
haviour of the elliptic flow in Au+Au collisions in the 
overall pt range is well reproduced by the model. In 
particular, a striking feature of the QGSM is that it is 
able to describe qualitatively the different px dependence 
of the mesonic and baryonic elliptic flow and reproduces 
a crossing of the elliptic flow for mesons and baryons at 
Pt > 1-4 GeV observed in the PHENIX experiment at 
RHIC. The centrality dependence of the integrated ellip- 
tic flow of charged hadrons in the QGSM agrees almost 
perfectly with the PHOBOS experimental data. This 
fact reflects the dominance of low px hadrons. 



In conclusion, the microscopic quark-gluon string cas- 
cade model based on the colour exchange mechanism 
for string formation is able to describe qualitatively and 
quantitatively well the bulk properties of the directed 
and elliptic flow measured in t^/snn = 200 GeV Au+Au 
collisions at RHIC. The limitations of the model should 
not permit a perfect agreement for all observables. The 
only signal which is really underestimated by the QGSM 
is the elliptic flow of hadrons with high transverse mo- 
menta pt > 1 GeV/c, although the t>2 of hadrons ob- 
tained within this model is already stronger than that 
of string models based on the FRITIOF scenario of ex- 
citation of longitudinal strings. The most plausible ex- 
planation of this discrepancy is an anisotropic character 
of jet absorption in a hot and dense asymmetric partonic 
medium. This effect, colloquially known as jet quenching, 
is intensively studied now Other processes related 
to the quark-hadron phase transition can be considered 
as well. The collective flow of hadrons appears to have 
a complex multi-component structure caused by rescat- 
tering and absorption processes in a spatially anisotropic 
medium. Therefore, further detailed investigations of the 
freeze-out scenario for hadrons are very important to un- 
derstand properly the flow formation and its evolution. 
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